Scale-dependent competing interactions: 
sign reversal of the average persistent current 
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The interaction-induced orbital magnetic response of a nanoscale system, modeled by the persis- 
tent current in a ring geometry, is evaluated for a system which is a superconductor in the bulk. 
The interplay of the renormalized Coulomb and Frohlich interactions is crucial. The diamagnetic 
response of the large superconductor may become paramagnetic when the finite-size-determined 
Thouless energy is larger than or on the order of the Debye energy. 
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Introduction. Rcnormalization, often accomplished us- 
ing the renornialization-group method, is one of the basic 
concepts in physics. It deals with the way various cou- 
pling constants (e.g., the electron charge or the Coulomb 
interaction) change as a function of the relevant scale 
for the given problem. The scale may be the resolution 
at which the system is examined, determined by its size 
and by the relevant energy for the process under consid- 
eration. Often, one knows the coupling constant's "bare 
value" at a less interesting (e.g., a very small or a very 
large) scale, and what is relevant for experiments is the 
value at a different, "physical" scale; the coupling on the 
latter scale is then used for the relevant physics. 

An important and particularly interesting application 
of this idea is in the theory of superconductivity. There, 
two competing interactions exist: the repulsive Coulomb 
interaction, starting on the large, microscopic energy 
scale-typically the Fermi energy/band width Ep, and the 
attractive phonon-induccd interaction, operative only be- 
low the much smaller Debye energy wd- By integrating 
over thin shells in momentum (or energy) space one 
finally gets the variation of any coupling g, being repul- 
sive or attractive, between a high-energy scale, a;>, and 
a low one, a;<, 
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(We use ?i = c = fee = 1.) Notice that 
pulsive/attractive interaction is "renormalized 
wards/upwards" with decreasing energy scale, 
makes superconductivity possible is that at lod the re- 
pulsion is much smaller than its value on the microscopic 
scale. Then, at wd the attraction may win and at lower 
energies the total interaction increases in absolute value, 
until it diverges at some small scale, which is just the con- 
ventional "Tc" of the given material. The upper panel of 
Fig. [1] presents this qualitative text-book picture. In 
this Letter we present quantitative results for a diffusive 
system whose size L exceeds the elastic mean-free path 
£ = v-pTi {v-p is the Fermi velocity and is the mean- free 
time). 



An intriguing situation offered by mesoscopic super- 
conductors, is the possibility to observe the "running", 
scalc-dcpcndent relevant interaction Q. Obviously the 
rcnormalization has to be stopped at the scale corre- 
sponding to the system size, which turns out for many 
physical observables to be the Thouless energy~the in- 
verse of the time it takes an electron to traverse the fi- 
nite system-which for diffusive systems is Ec = D/L^ 
[D = wp^/S is the diffusion coefficient) Moreover, 
Ec is also the energy scale for many mesoscopic phenom- 
ena. Thus, by measuring a quantity which depends on 
the interaction, as a function of L, one may deduce the 
scale-dependence ("running" value) of that interaction. 
When that quantity is sensitive to the sign of the inter- 
action (e.g., a linear function, up to corrections, of the in- 
teraction, often resulting from perturbation theory), one 
may observe a rather dramatic change of the sign of the 
interaction-dependent term in that quantity! 

A unique example for the latter is the orbital magnetic 
susceptibility of a mesoscopic (or even nanometric) sys- 
tem. Most convenient for calculations is the persistent 
current Q flowing in a mesoscopic normal-metal ring in 
response to an Aharonov-Bohm flux through its opening. 
The magnetic moment induced by very small fluxes due 
to these currents is expected to be of the same order of 
magnitude as that of a small disc made of the same mate- 
rial, and having similar radius and thickness, in response 
to the same flux. 

Moreover, it is well known 0, [H| that the impurity 
ensemble-averaged persistent current for noninteracting 
electrons is smaller by typically two orders of magnitude 
than the interaction-dominated terms. Therefore the lat- 
ter represents rather well the whole magnetic response of 
the electrons and its measurement as a function of L 
will amount to monitoring the renormalized interaction. 
This picture is borne out by our calculations. We imple- 
ment the ideas of Refs. Handi in the derivation of the 
superconductive fluctuation-induced partition function, 
and obtain the effect of the two competing interactions 
in the Cooper channel. We use this result to study the 
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average persistent current of a large ensemble of metallic 
rings in the presence of these interactions. 
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FIG. 1: Upper panel: the renormalized interactions, Eq. 
([T]), as functions of the energy, for a metal which is a su- 
perconductor in the hulk. When only the attraction exists, 
its absolute value is renormalized up (with aj> = ljd) un- 
til it diverges at T° (the dash-dotted curve). The repulsion 
(dashed curve) is renormalized down from its bare value at 
cj> = Ey- The total interaction starts at uiu from the bare 
value ga{oJY)) + gr(oju) and then its absolute value is renormal- 
ized up until it diverges at the "true" transition temperature 
Tc (full curve). In a mesoscopic system the renormalization 
stops at the Thouless energy and therefore the relative loca- 
tion of Ec vs. oiD and the ratio |(ja|/gr(max{i?c,<^D}) will 
determine the sign of the effective interaction at ojd- The fig- 
ure is drawn for demonstration purposes with the unrealistic 
parameter Ef/luo — 2.5. Lower panel: the lines in the plane 
Qr — Ec/cdD where the first harmonic of the current changes 
sign from paramagnetic (above the curve) to diamagnetic (be- 
low), for three values of the attraction Qa. Here and in all 
figures below, T = a;D/125 ~ 2 - 4K, and Ef/ujj^ = 25. 

The renormalization of the interactions has to be 
carried out down to an energy scale roughly given by 
max{T, Ec}, where T is the temperature. We denote 
by ga the attraction and by gr = [g^^ + ln(£'p/wp)]^^ 
the repulsion, both at wd, such that for \ga\ > gr, the 
material in its bulk form is a superconductor. When 
Ec > T then depending on the ratios Ec/ujd and 
|(7a|/5r(max{i^c, wd}) the resultant interaction may be 
repulsive or attractive, and the response may be para- 
magnetic or diamagnetic. These two types of behavior 
are reflected by the sign of the average current at small 
fluxes. The lower panel of Fig. [T] presents a "phase di- 



agram" showing the regions where the first harmonic of 
the current is paramagnetic or diamagnetic, as a function 
of Ec/oj-£,. Note in particular the existence of a region 
where although the attraction \ga\ exceeds the repulsion 
gr, the current is paramagnetic, contrary to the usual ex- 
pectation 0. The Thouless energy can be increased by 
making the ring smaller, and/or less disordered. 
Renormalizing the two interactions. The Hamil- 
tonian we consider consists of a single-particle part, 
i/o(r) = (— «V — eA)^/(2m) — ^ + u{r), and a term de- 
scribing local repulsion and attraction, of coupling con- 
stants gr = XrM/V and ga = XaAf/V, respectively, 

+ (A,-|AJ)^|(r)^j(r)^^(r)^^(r)) , (2) 

where N' is the density of states, V the system volume, 
and tjj^ {tpt) destroys (creates) an electron of spin com- 
ponent a at r. Here, Ijr is the bare repulsion, /i is the 
chemical potential, u{r) is the disorder potential due to 
nonmagnetic impurities, and A is the vector potential, 
with A = ^/ L, ^ giving rise to a persistent current. As 
the magnetic fields needed to produce such currents are 
rather small, the Zeeman interaction will be ignored. 

The partition function, Z, corresponding to the Hamil- 
tonian ^ is calculated by the method of Feynman 
path integrals, combined with the Grassmann algebra 
of many-body fermionic coherent states Introducing 
the bosonic fields </'(r, r) and A(r,T) via two Hubbard- 
Stratonovich transformations one is able to integrate over 
the fermionic variables; the result is then expanded to 
second order in the bosonic fields M 
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with the action S given by 

S = Y,[f3M\A^jy\gJ+(im<i.AV9r 

- (Aq,. + *0q, J (a;, + J Hq,,) , (4) 



and Zq denoting the partition function of noninteract- 
ing electrons (/3 = 1/T). Equation ^ represents the 
partition functioii^ due to superconductive fluctuations, 
with the action S, a sum over the wavevectors q and 
the bosonic Matsubara frequencies v = 2TrnT, written 
in terms of the polarization 11 averaged over the impu- 
rity disorder It is this object in which the different 
renormalization of the competing interactions manifests 
itself. 

In a diffusive system, the disorder-averaged polariza- 
tion is given as a sum over fermionic Matsubara frcquen- 
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cies, Lo = 7rT(2n+ 1), 
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where 8 is the Heaviside function. The first of Eqs. ([5]) 
is vahd for |2a; + i/j <^ and qi <^ 1. It represents a 
divergent sum and therefore a cutoff is needed. The cutoff 
energy is determined by the relevant interaction. The 
sums in the terms of Eq. ^ involving A are bounded 
by uju and therefore are cut off at /(27rT) - 

1; these sums are denoted 11"°. The sum multiplying 
is cut off by Ep, i.e., n„iax — Ep/{2'!TT) - 1; this 
sum is denoted by . Therefore, there appear two 
polarizations, the one cut off by and the other by E-p 
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where ^ is the digamma function. 

Performing the integration over the fields yields 
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where we have introduced the wavevector- and frequency- 
dependent renormalized repulsion ^^(q, J^) 
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In the limit :s> {Dq'^,T} one may neglect F^^^, in the 
arguments of the digamma functions and retrieve (ex- 
ploiting the asymptotic form of VP) the well-known result 
= + ln(iJp/wj-,). In a finite ring, however, the 
smallest wavevector is of the order 1/L and therefore at 
low temperatures this procedure is valid only when the 
Thouless energy is smaller than the Debye energy, i.e., 
for large enough (and/or clean enough) rings. Carrying 
out the integrations over the fields A in Eq. (O yields 
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Equation ([9]) is our result for the superconductive 
fiuctuation-induced partition function. It is based on the 
renormalized repulsion, Eq. 

The fiuctuation-induced average persistent cur- 
rent. The thermodynamic persistent current is ob- 
tained by differentiating the free energy with respect 
to the dimensionless flux, Lp = $/<I>Oi i-e., I = 
[eT) / {2Ti)d{hi Z) / dip where $o = /i/e is the fiux quan- 
tum. In the ring geometry, the flux enters the longi- 
tudinal component, q\\, of the wavevector q, q^^^ = q^ = 
2'!T{n + 2ip)/ L where n is an integer Hence (exploiting 
the Poisson summation formula) 



/ = sin(47rrn(p)/„ 



(11) 



where 
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dxe'-"dHa,^^^,]/dx , (12) 



with ax,m,iy [Eq. (fTOj) ] given in terms of Fx^m.v ~ 
0.5 + \v\/{AttT) + {TTE^x'^)/{m^T) [see Eq. ©]. The 
integral in Eq. (|12p is handled by a contour integration, 
and the sign of the result which determines the magnetic 
response, depends on a subtle balance between the con- 
tribution of the poles and that of the zeros Q . 

It is customary in studies of persistent currents to con- 
centrate on the lowest harmonics, notably on /m=i, since 
higher harmonics are exponentially smaller compared to 
the lowest ones and are also more sensitive to inelastic 
scattering (dephasing) and pair-breakers (e.g., magnetic 
impurities) Q. However, when T < Ec~ values are con- 
sidered the relative magnitude of harmonics higher than 
the first one increases, and their contribution should not 
be discarded. 

When only one of the two competing interactions is 
accounted for, the nature of the magnetic response is dic- 
tated by the sign of that interaction, being diamagnetic 
for attraction and paramagnetic for repulsion When 
both interactions are present the picture is changed, as 
is exhibited in the lower panel of Fig. [U obtained by 
numerically calculating Im=i as a function of the Thou- 
less energy Ec- The sign of Im=i depends on the ratio 
Ec/uj-D- As the latter increases, it may change from nega- 
tive to positive. On the other hand, when g^ > \g^\ (and 
the bulk material will not become a superconductor at 
a finite temperature) the current remains paramagnetic 
for all values of Ec/oj-d- 

As mentioned, it is when T > Ec that the first har- 
monic represents faithfully the persistent current. Figure 
[2]depicts the first three harmonics of the current, as func- 
tions of the Thouless energy (both normalized by ujd)- 
The figure indicates that the value of Ec at which the the 
sign of Im is reversed increases with m. Figure [3] demon- 
strates the dependence on the strength of the repulsion. 
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FIG. 2: Im{Ec) for m — 1, 2, and 3 (see the legend), with 
= 0.2 and 3, = 0.1. 

As it increases the Ec/uju^"cangc in which the response 
is diamagnetic diminishes, and so does also the maximal 
value of the diamagnetic current. In the paramagnetic- 
response regime the amplitude of the current increases 
with Ec Q • Figure H] shows the flux dependence of the 
current, for two distinct values of gr, and displays the 
possibility for a sign reversal of the magnetic response at 
low flux. One notes in Fig. [2]that at Ec/iou = 0.76 /m=i 
is small compared to /m=2,3; the dashed curve in Fig. |4] 
reflects the dominance of the higher harmonics over the 
first one for relatively small gr- For a higher value (solid 
line there) the first harmonic is dominant. 



traction suffices to lead to a superconducting phase in 
the bulk material, we find that it does not ensure a dia- 
magnetic response of the mesoscopic system. The reason 
being that in mesoscopic rings the orbital magnetic re- 
sponse owes its very existence to a finite Thouless energy. 
Therefore the latter, when large enough can cause the 
response to be paramagnetic, albeit the strength of the 
attractive interaction. As the Thouless energy may be 
controlled experimentally, one may hope that the predic- 
tion made in this paper will be put to an experimental 
test. Relatively large values of persistent currents may 
be achieved in molecular systems, and small discs are ex- 
pected to behave similarly at small fiuxes. Finally we 
remark that there may already be experimental indica- 
tions to the validity of our prediction. Reich et al. 
found that thin enough gold films are paramagnetic. This 
may well be due to a small grain structure. 
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FIG. 3: I{Ec) for three values of the repulsion and \ga\ = 0.2. 
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FIG. 4: The current as a function of ifi. Solid (dashed) curve 
is for = 0.17 {fjr = 0.1), with \ga\ = 0.2 and Ec/uju = 0.76. 

Discussion. We have obtained the wavevector- and 
frequency-dependent renormalizcd repulsion, Eq. ([5]). 
The two cutoffs on the bare attraction and repulsion, 
cjq and Ep, respectively, conspire to determine the sign 
of the magnetic response. Whereas an overall total at- 
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